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1. INTRODUCTION 

Numerous types of electric generators such as the linear, induction, permanent magnet and 
synchronous generators are used together in order to achieve the electric power generation. The latter ones 
are utilized almost exclusively in power systems as a source of electrical energy due to its accessibility for 
regulation. 

Synchronous generator constant voltage at the generator terminals is necessary for satisfactory main 
power supply [1]. The output voltage can be affected by diverse disturbing factors such as random 
fluctuations in load, fault conditions. For this purpose, automatic voltage regulators (AVRs) are utilized for 
guarantying a stable terminal voltage of the synchronous generators by adjusting the excitation voltage. 

Unfortunately, without any controller, the AVR system does not provide a suitable performance [2]. 
For this, several control methods have been studied and developed in the literature [3-5]. One of them is the 
proportional-integral-derivative (PID) and (PI) controllers which are widely used in industrial control 
applications and science like AVR systems [4], [6]. This popularity is due to their effectiveness, ease of 
implementation and their simple structure [7]. 
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However, the tuning aspect of PID parameters is not easy for researchers and plant operators [8]. In 
this context, many different methods have been suggested to determine the parameters of PID controller. The 
most well-known method used in this target is the Ziegler-Nichols approach [9, 10]. Nevertheless, this 
method is complex and time consuming for systems with non-linearities such as synchronous generator [11]. 
Generally, it delivers a huge overshoot in the response of AVR system and is not suitable to achieve a fast 
dynamic response [10], [12]. 

In order to remedy these problems, different optimization methods have been adopted for tuning the 
optimal parameters of PID controller of an AVR system. Such algorithms included genetic algorithm (GA) 
[13], particle swarm optimization (PSO) [4], invasive weed optimization (IWO) algorithm [14], Artificial 
Bee Colony algorithm (ABC) [6], cuckoo search (CS) [10], and others. 

Furthermore, an effective control of the terminal voltage of the synchronous generator requires an 
accurate model. Most of researchers have used transfer functions to represent the model of AVR system 
when the generator model is also represented with a transfer function, which differs between them to the 
order of the model. A low order (first and second-order) is used in [4], [13] and a high order is adopted in the 
works [15], [16]. To determine the transfer functions, a linearization 1s needed. In other words, a limited 
development is performed around the operating point. Therefore, the model is validated only around the 
operating point. 

In this work, for output voltage regulation of synchronous generator, the nonlinear model of the SG 
is used directly without needing a linearized mathematical model of the generator. Also, a varied reference 
voltage is adopted in order to decrease the overshoot resulting in the response of the terminal voltage. For 
this study, a PI controller is used to vary the duty cycle of the DC chopper of step down type for controlling 
the excitation voltage and therefore the terminal voltage of the generator. In addition, the PSO algorithm is 
applied for obtaining the parameters of the PI-Controller. The performance of the PSO-PI controller is tested 
with the GA-PI controller. To verify the validity of the varied reference voltage, a comparison of the 
response of the terminal voltage using the varied and the fixed reference voltage is carried out in this paper. 


2. MODELLING OF SYNCHRONOUS GENERATOR 

The study was carried out on a salient-pole synchronous generator. Figure | shows the simplified 
structure of a salient-pole synchronous generator, which comprises three armature windings, one excitation 
or field winding and short circuited damper winding in each axis (d-axis and q-axis). In order to obtain the 
fundamental equations of a synchronous machine, the principal assumptions adopted are [17]: 
- The three-phase stator windings are symmetrical and are sinusoidally distributed in space around the 

air-gap. 

- Hysteresis effects are ignored. 
- There is no magnetic saturation. 
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Figure 1. Synchronous generator windings with dampers. 


Considering the generator convention for the stator windings, in the Park framework, the dynamic 
equations of a synchronous generator related to a reference linked to the rotor are written, as follow [18]: 
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d d d 
Va = —Rsla + ava — WW; Vg a —Rslg + q Pa + wpa; Vr = Rplp + a PF 
d d 
In which;V,R,/J,@, and p denote voltage, resistance, current, angular speed and flux 
linkage.s, F, D, Q indicate armature, field and damper windings. Subscripts d and q stand for the direct and 
quadrature axes, respectively. The Flux linkage relationships are expressed as (2): 

Pa = —Lala + Mralr + Mpalp; Pa = —Lqlq + Moqalgi We = Lele + Meplp — Mrala; 

Wp = Lplp + Marple r= Moala; Wo = Lolo z Moqlq (2) 
where; Lp, Lp and Lg are the inductances of the main field, direct and quadrature damper windings 
respectively. Lg and L, are the inductances of the d-axis and q-axis stator winding. Mpa, Mra, Mpa and 
Mogare the mutual inductances between the stator d-axis and the rotor main field, the rotor main field and the 
direct damper, the d-axis stator and the direct damper, finally between the q-axis stator and the transverse 


damper respectively. 
The dynamic equations of the generator are given as follow (3): 


d d d 
Va = =R;l4 + WLglg T WMoglo = La id + mea ae + gee 
Va = —@Lglg — Rslq + OMepgle + OMpglp — Lg ala + Moq zle 
d d d 
Vr = Relr a Mra aid + Lr Ti |F + Mpp PrE 
d d d 
0 = Rplp + “he ae + a E + Lp aD 
2.1. Generator connected to an inductive load 


Let’s consider that an inductive load is connected to the generator, after Park’s transformation, the 
direct and quadrature components of the load voltage are represented by the following (4): 


Va] [—-R, oly] fla -L O 14 fla 
vi 7 Bex a | cr | 0 e py (4) 


where; R; and L, are the resistance and the inductance of the load. 
By combining the (3) and (4), the global equation of the generator is given by (5): 


[V] = [RIU] + [L]—— (5) 


where, 


[]= [a l4 Ir lb lę] ,[V]=[0 0 Vp 0 0J 


RL-R, œ(La-L) 0 0  -wWMggq Ly — La 0 Mra Mpa 0 
o(L,—Lg)  R—R,  @©Mpgqg ©Mpa 0 0 Ly —Lg 0 0 Mag 
[R] = 0 0 Ry 0 0 ; [L] =| —Mra 0 Le Mpp 0 
0 0 0 RS 0 —Mpa 0 Mrp Lp 0 
0 0 0 0 Ro 0 —Magq 0 0 Lo 


So the state Equation is (6): 
-7 = [L] [R] U] + [L] >V] (6) 


The stator terminal voltage V,,(t), Vps (t) and V.,(t) can be derived numerically from the values of V4 (t) and 
V,(t) by using the inverse Park transformation defined by (7) according to the machine’s electrical angle 
0 [17]: 
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cos 0 — sin 0 = 
P(@)"' = = 2] cos cos (o — =) — sin (o — =) = (7) 
cos ( =) — sin (o + =) = 


Such that Vape = PO) * Viass by considering that the output voltage is well balanced, the homopolar 
component is null and by tacking, 0 = wt + 6, the stator terminal voltage expressions V,,(t), V,,(t) and 
V.,(t)are given by: 


V,,(t) = iE (Va(t) cos(wt + 69) — V,(t) sin(wt + 69)) (8) 
V, (t) = i (Va(t) cos (wt + 0) —) - VC sin (wt + 6 - =)) (9) 
V- (t) = i (va (t) cos(wt + 0o + =) — V(t) sin(wt + 0o + m) (10) 


where, 
w : machine rotor speed in electrical radians/seconds. 
0o: is the angle between the stator of phase ‘a’ and the rotor field axis (d axis) at the instant t=0. 


3. OUTPUT VOLTAGE REGULATION OF SYNCHRONOUS GENERATOR 
Figure 2 describes a block diagram to ensure the output voltage of the synchronous generator at the 
required value during load variations. 


Reference voltage Vef (t) 














Figure 2. Block diagram of the proposed automatic voltage regulator (AVR) with PI-controller 


This block diagram contains a controller, a converter, a main power source of the converter and a 
synchronous generator. The generator field winding requires DC supply and it is adjusted by a DC chopper 
(buck type). The terminal voltage of the generator (V,(t)) is compared with the reference voltage (V,ef(t)) 


to give the error (e(t)) as shown in (11). 
e(t) = Vref t) — Vie) (11) 


The duty cycle a of the DC-chopper is regulated by the PI-controller and the output voltage V; of the 
DC-chopper is given by(12) depending on the duty cycle æ [19]. 


Vy = a. Vinput (12) 
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Where; 
V; is the output of the DC-chopper or field voltage (ve) and Vinputis the supply voltage to the DC-Chopper. 
By controlling the duty cycle a, the field voltage is adjusted to regulate the output of the generator. 


3.1. Proportional Integral controller (PI Controller) 
The PI controller consists of two parameters, the proportional gain, denoted K, and the integral gain, 
denoted K;. Figure 3 shows the typical structure of PI controller. 


















r(t) _ e(t) = 
+ 
7 t 


Figure 3. Basic block of proportional integral controller (PI controller) 


The PI controller in the time domain is described by the relation [20]: 
u(t) = Kye(t) + K; f, e(t)dt (13) 


where, e(t) and u(t) are tracking error and the control signals represented in the time domain. 
The transfer function of this controller is represented in the following (14): 


u(s) = (k, + < e(s) (14) 


4. TUNING OF PI-CONTROLLER FOR OUTPUT VOLAGE REGULATION USING PSO 
APPROACH 

4.1. PSO Algorithm 

Particle swarm optimization (PSO), developed by Eberhart and Kennedy [21] in 1995, is aswarm 
intelligence meta-heuristic inspired from the group behavior of animals, for example movement of bird 
flocks or fish schools [22] can be simulated. The algorithm starts by initializing a population of random 
solutions called particles and searches foroptima by updating generations through the following velocity and 
position update equations. 
The new velocity of each particle is formulated by the relation: 


vi(k +1) = wv,(k) + cr; x (pbest; (k) — x; (k)) + CGT, X (gbest(k) — Xi (k)) (15) 
The position of each particle is updated in each generation according to the following expression: 


In (15) and (16), v;(k + 1) is the new value of velocity vector, v;(k)is the previous velocity vector, 

x;(k + 1)is the new position vector, x;(k) is the previous position vector, pbest;(k) is the particle's local 
best known position, and gbest(k) is global best known position. c, and c, are positive constants called 
respectively coefficient of the self-recognition component and coefficient of the social component. w is the 
inertia factor, 7, and r, are random numbers between zero and one. 

The parameter w shall be properly set such that the impact of the previous velocity and the 
magnitude of the position shifts depends on the values of the parameters c,and c3. Therefore, the choice of 
these parameters is very important for the convergence of this method. In the original PSO algorithm, the 
positive acceleration constants c,and c, are usually set between 0.5 and 2 [21]. 


4.2. Implementation Of A PSO-PI controller 


The block diagram of PI controller tuned by PSO algorithm is shown in Figure 4. The principal 
objective of the PSO algorithm is to search the optimal value of K, and K; based on the objective function F, 
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which is given by the relation (17). In this article, F represents the quadratic output error between the 
reference voltage and the terminal voltage calculated from the adopted model. 


F = eV — V)” (17) 


Where n is the number of sampled points, V,efiis the reference voltage, and Vis the terminal voltage of the 
generator. 





Process 


Figure 4. Block diagram of PI controller tuned by PSO 


where; 
Vinput : is the supply voltage to the DC-Chopper. 
Ver =220v (r.m.s) is the fixed reference voltage. 

In this paper, the terminal voltage is the max of the stator terminal voltage V,,(t). Note that the latter 
one is obtained by a numerical resolution of the non-linear system (6) by means of the fourth-order Runge- 
Kutta method using the parameters of the reference [23]. Figure 5 shows the flowchart of PI-PSO approach to 
regulate the terminal voltage of a synchronous generator. The parameters for the PSO algorithm are collected 
in Table 1. 


Initialize PSO parameters randomly 
Call PI Controller 
Evaluate objective function for initial population by (17) 


Update velocity and swarm position (move ag 1 
the particles with the new PSO velocity) ) Iter =Iter-+ 1 | = Iter + 


Call PI Controller No 
Fitness<Tolerance 
Evaluate objective function for eachparticle by (17) 


Yes 
Compute pbest, gbest 
Calculate new searching direction 


Get best solution for PI parameters (K, and K;) 





Figure 5. Flowchart of the PSO-PI controller 


Tablel. Parameters of PSO algorithm 





PSO propert Value 
Number of dimensions of the search space 2 
Initial particles positions Random 
Population size 10 
Inertia weight factor w 0.8 
Acceleration constant c4 1.6 
Acceleration constant c, 1.5 
Tolerance 1076 
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5. SIMULATION RESULTS AND DISCUSSION 
Three different loads are applied to the generator at different instant in order to show the impact of a 
load on the terminal voltage. Also, the generator loading conditions are: 
- The generator starts with the inductive load L, characterized by R,, = 3002; L,, = 0.1H 
- After one (01) second, the load L,(R,2 = 10.2; L,, = 0.2H) takes the relay. 
- After (02 seconds), the load L3(R,3 = 80.2; L;3 = 0.1H) is applied. 


5.1. Before regulation (without PI-controller) 

In this case, the field excitation voltage must be kept constant and the generator operates at rated synchronous 
speed (w = 157rad/s). Figure 6 shows the terminal voltage response of the SG before regulation (without 
PlI-controller) for different loading conditions. 

















oO 0.5 1 1.5 2 2.5 3 3.5 
Time(s) 


Figure 6. Output voltage response of the SG before regulation for different loading conditions. 


As shown in Figure 6, the variation of the load connected to the generator has a direct influence on 
the voltage generated. Furthermore, the response time is large, as well as the voltage drop and the 
disturbances caused by the load changes. 


5.2. After regulation 

In this part, the reference voltage is taken at the value (Vref=220 v (r.m.s)). The regulation of the 
output voltage of the synchronous generator using the (PSO) algorithm is simulated by interaction between 
MATLAB and Simulink. In order to evaluate the performance of the PSO-PI controller, the same simulation 
is implemented using the genetic algorithm (GA) [24]. The parameters used in this study for the GA-PI 
controller are summarized in Table2. Note that, the lower and upper bound of the gains (Kp and Ki) are the 
same for both algorithms. Table 3 provides the optimal gain parameters of the PI-controller tuned by GA and 
PSO algorithms. Figure 7 shows the superposition of the terminal voltage response of the SG with PI 
controller tuned by GA and by PSO through diverse loading conditions. The parameters (Kp and Kj) 
evolution versus the number of iterations is illustrated in Figure 8. 


Table 2. Parameters of GA algorithm 


GA property Value 
Number of dimensions of the search space 2 
Initial population Random 
Population size 10 
Crossover probability 0.75 
Mutation probability 0.06 
Tolerance 1076 


Table 3. PI controller gain values. 


Controller type Kp K; Number of iterations Time computing (sec) 
PSO-PI 0.00815 0.72149 186 5252.8656 
GA-PI 0.00809 0.71850 313 11428.5708 
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Figure 7. Superposition of the terminal voltage response of the SG with PI controller tuned by GA and by 
PSO through diverse loading conditions. 
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Figure 8. Parameters evolution of the PI-controller: (a) Kp, (b) Ki 
(solid red line: PSO, solid blue line: GA) 


Output voltage regulation of synchronous generator using PSO algorithm-based PI controlle (Fatiha Habbi) 


1224 O ISSN: 2088-8694 


As it can be seen in Table 3, the proximity of the PlI-controller parameters (Kp and K1) confirms 
well the convergence of the (GA) and the (PSO) algorithms. Furthermore, the PSO-PI controller consumes 
less time and has less iterations than the GA-PI controller as shown in Figure 8.However, the Figure 7proved 
that the overshoot resulting in the response of the terminal voltage of the generator is bigger. 

In this study, The PSO algorithm is used to determine optimal parameters (Kp and Ki) of the PI- 
controller, which regulates the duty cycle of the DC-chopper, as shown in Figure 9 to ensure the convergence 
of the terminal voltage of the generator to the desired value. The generator armature current in phase (a) 
during different loading conditions is plotted in Figure 10. 
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Figure 9.The duty cycle of the DC-chopper Figure 10. Generator 
armature current (Ias (t)) 


5.3. Imrovement of the overshoot 

In order to reduce the overshoot resulting in the terminal voltage response of the SG, we adopted a 
varied reference voltage. In other words, the previous reference voltage (Vref=220 v (r.m.s)) is replaced with 
a ramp form during tọand continued with a constant voltage during t — tọ as shown in Figure 11. This 
improvement takes into consideration that the reference voltage does not suddenly go from zero to the 
maximum value of the set point, but it increases progressively until the voltage set point is reached. 
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Figure 11. The varied reference voltage 


where; 
to: is the ramp time corresponding to the time required to reach the generator voltage set point. 
t, : is the simulated time. 

The varied reference voltage can be expressed as (18): 


220 
Vier = (P)er<t, (18) 
220ty <t<t, 

In this case, tọ takes the value 0.06s. After giving the same parameters settings above to the PSO 
algorithm, optimal PI controller parameters were obtained as shown in Table4. The comparative output 
voltage response of the generator using the PSO-PI controller with the fixed and the varied reference voltage 
is Shown in Figure 12. 
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Figure 12. Comparison of the terminal voltage response 


From the comparison shown in Figure 12, it can be seen that the terminal voltage response with 
varied reference voltage has small overshoot and short settling time compared to the terminal voltage 
response with a fixed reference voltage. The transient response specifications (overshoot percentage, the 
settling time and the rise time) are calculated using the relationship provided by the reference [25]. The 
comparative results obtained are summarized in Table 5. 


Table 5. Comparison between the output voltage of the synchronous generator using 
the varied and the fixed reference voltage 


Parameters Varied reference voltage Fixed reference voltage 
Overshoot percentage % 1.89 6.59 
Settling time (s) 0.29 0.34 
Rise time (s) 0.02 0.02 
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6. CONCLUSION 

This work presents a regulation of the output voltage of a synchronous generator using the nonlinear 
model of the SG. The PSO algorithm is utilized for the determination of the optimal parameters of the 
proportional integral (PI) controller. The latter one is used to adjust the duty cycle of the DC chopper for 
controlling the field voltage and consequently the output voltage of the generator. The performance of the 
PSO-PI controller is tested with the GA-PI controller. In terms of computational time and iteration number, 
the obtained results show the superiority of the PSO-PI controller versus the GA-PI controller but both the 
PSO-PI controller and the GA-PI controller produces a bigger overshoot in the terminal voltage response of 
the generator. To circumvent this problem, a varied reference voltage is chosen. The comparison between the 
terminal voltage response using the fixed and the varied reference voltage proves that the latter one has the 
capability to reduce the overshoot and also the settling time. 
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